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Abstract 

Coupled plasmon modes are studied in graded plasmonic waveguides, which are periodic chains of metal- 
lic nanoparticles embedded in a host with gradually varying refractive indices. We identify three types of 
localized modes called "light", "heavy", and "light-heavy" plasmonic gradons outside the passband, ac- 
cording to various localization. We also demonstrate different transitions among extended and localized 
modes when the interparticle separation d is smaller than a critical d c , whereas the three types of localized 
modes occur for d > d c , with no extended modes. The transitions can be explained with phase diagrams 
constructed for the lossless metallic systems. 
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Recent explorations of surface plasmon (SP) waves sustained by metallic nanostructures have 
been successful in nanooptics, see, e.g., Refs. 1-2. Specifically, guiding and manipulating electro- 
magnetic (EM) energy below the diffraction limit with SPs have attracted much effort, concentrat- 
ing on routing plasmon signals in periodic noble metal nanoparticle chains . 3 i 4 i 5 ' 6 In these so called 
plasmoic waveguides (PWs), metal nanoparticles can sustain resonant collective oscillations (plas- 
mons) of their conduction electrons. When driven by an external light field, these electron oscilla- 
tions couple to the visible optical excitation to form SPs bound around the nanoparticles, therefore 
possess a lot of novel properties . 1 i 2 i 3 i 4 i 5 i 6 The dispersion relations of coupled SP modes in these 
PWs are like those of propagative EM waves in one dimensional Bragg gratings, 4 rather more like 
electron dispersions within "tight-binding" atoms in a solid. The plasmon resonant bands, which 
are quasicontinuous for finite system and continuous for infinite system, are determined by the 
interparticle couplings and almost centered by individual particle's site Mie resonances^ which 
are, crucially dependent on a rich set of parameters such as the nanoparticle's size, shape, as well 
as surrounding medium.^ 

In PWs, however, longitudinal confinements of the plasmon excitation should also be of great 
interest and have been explored little so far. Here we propose to use a host with graded refrac- 
tive indices along the nanoparticle chain — graded plasmonic waveguide (GPW). This belongs to 
a subclass of deterministic nonperiodic media which has recently been revealed to exhibit inter- 
esting spectrum transitions in both elastic and plasmonic systems. 8 - 9 As compared to the case of 
homogeneous host (e.g., vacuum), apparently a graded host has two consequences: (i) shift of 
Mie resonance of individual particle, simply because dielectrically denser host results in a heavier 
effective optical electron mass m* in the metallic particle;- (ii) screening effect on the interparticle 
couplings. It turns out that GPWs are far more intriguing than expected; very rich longitudinal 
(axial) localization effects are discovered. 

We consider a linear chain of iV spherical metallic nanoparticles immersed in a dielectric host. 
We take lossy Drude form dielectric function e(u) = 1 — Up/u(uj + iT) for the nanoparticles while 
assume that the host has its dielectric constant varying from the chain's left-hand side to the right- 
hand side along x axis as e2(x n ) = ei + cx n /l, (n = 1, 2, • • • , N), where x n denotes the position 
of the n-th nanoparticles, / the total length of the chain, and c the coefficient of dielectric gradient. 
A point dipole model is believed sufficient in capturing the fundamental dispersion characteristics 
of coupled plasmons in PWs, 4 ? we expect this to be also true in GPWs. The coupled equation for 
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point dipoles reads 

G(w)p = 0, (1) 

where p is the Arrowed column vector of dipole moment oscillating with the frequency cu, 
G nn (uj) = l/p n (u),(n = 1,...,N) and G nm {uS) = -f nm , (m ^ n). Here f3 n (u) = 
^2(x n )[e(uj) — e 2 (x„)]/[e(o;) + 2e 2 (x n )] represents the polarizability of the n-th particle, and T nm 
denote the EM interaction on the n-th particle due to the m-th particle. 5 Equation O can also be 
regarded as an eigenequation with normal mode frequency uj a determined by det{G(u)} = 0, 
while to a generally being complex valued. However, it is hard to determine uj a by directly solving 
this equation which is a nonlinear complex transcendental equation with respect to cu. Alterna- 
tively, we employ the inhomogeneous equation G(uj)p = v to study the plasmon excitations, as 
inRef. 5. 

FigureUJshows, for the longitudinal polarization (along the chain axis x), the real part of dipole 
moments obtained via p = G~ 1 (uj)v, where G -1 ^) is evaluated at a real valued frequency uj 
indicated in every panel. 5 In these calculations, we have set all but the row of N/2 in v (a column 
vector) equal to zero. The GPW parameters are eg = 3.0, c = 1.0, A" = 100, and interparticle dis- 
tance do = 4a (3a) for the left (right) panels, where a is the radii of the nanoparticles. To highlight 
the graded effect, only near- field couplings (f nm oc 1/ \x n — x m | 3 ) are taken into account. We note 
that the fully retarded couplings were shown to be remarkably important in the case of transverse 
polarization. 5 6 However, their role in modifying the dispersion is a subject of controversy. 10 These 
issues become more ambiguous in the presence of loss (radiative or nonradiative), which needs to 
be addressed appropriately further. Here we mainly focus on the longitudinal polarization. 11 The 
demonstrated localizations shall hold even with full couplings included,^*^ in view of a recent 
exploration in the localized vibrations of optically (long-range) bounded structures. 12 It is seen in 
Fig. d that there exist several different solutions which correspond to modes excited by external 
periodic "forces" subject to the central site n = N/2 with different frequencies: (1) at relatively 
low frequency, modes are excited in the right-hand side (larger e 2 ), which implies the externally 
injected energy moves to the right [Figs. [Ha) and[TJd)]; (2) at intermediate frequency, it propagates 
to both sides [Figs. [Tib) andUJe)]; (3) while at high frequency, it tends going to the left [Figs.[JJc) 
and[TJf)]. We can regard these excitation modes as being localized or extended depending on their 
spatial extent. However, if T is large, we can not simply regard the solutions (not shown) in this 
way. This is because it is hard to distinguish the localization effect from the loss which leads to 
attenuation of the excitation in the medium.— 
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In order to gain more insights into the peculiar and principal characteristics of GPWs, we 
therefore focus on weak loss or lossless regime hereafter. It is possible to linearize Eq. © with 
respect to uo 2 

(F - u 2 I)p = 0, (2) 

where i" is identity matrix and F nm = u 2 — G nm (ui 2 ) / [ G' nn (uj^) . 9 Here the prime indicates deriva- 
tive with respect to to 2 and uj n is the dipolar Mie resonant frequency of the n-ih. particle. Equa- 
tion © represents a good approximation of Eq. © at u ~ cu n . In considering the case of weak 
loss, we can write the eigenspectrum {cu a = uj' a + icu'^,p a } and v = J2 a VaPa- If me frequency u 
is close to a eigenfrequency u> ao then significant contribution to p comes from a few eigenmodes 
whose eigenfrequency is close to cu. Let us call the set of modes S that has significant contribution: 
\u>' a — uj' ao \ <C oj'u for a G S. Thus, the real part of the response function (e.g., plotted in Fig.Q) 
is Re[p] = Y.aes^iuj ~ u}' a )p a /2u}[(u - uj' a ) 2 + This means by solving p = G _1 (u;)v 
for weak loss in the Drude model, the solutions we get in Figs . E a )-Ef) are ver Y close to a single 
complex eigenmodes in the GPWs. 

In the linearized lossless case, the coupled plasmon dispersion can be directly obtained by 
diagonalizing F, resulting in real eigenspectrum {u a ,p a }. 9 The results are shown in Figs. Eta) 
andEtd) (solid lines). Still these are in the near-field coupling approximation, for the same GPW 
structures studied in Fig. [T] except for T = now. Also the density of states (DOS) D(uS) = 
^ a 5{uj — uj a ) /N are shown in Figs.|2tb) and|2e). The real valued eigenmodes p a (not shown) 
have quite different spatial extents, much like the excitation patterns in Fig.[l] This is also reflected 



by their inverse participation ratio 9 (IPR) ^ n=1 pt( n ) I E n =i Pa( n ) shown in Figs. |2c) and 
Etf). The IPR is of the order of 1/N if modes are extended over the system with N degrees of 
freedom, while it becomes larger than 1/N for spatially localized modes. We see that relatively low 
frequency modes are localized at the side of larger e 2 while high frequency modes are localized at 
the side of smaller e 2 . We name these modes as "heavy gradons" and "light gradons", respectively. 
More interestingly, coupled plasmon modes at intermediate frequencies between two transition 
points cu^do, c) < uj < o;h(c?o, c) are extended 9 when d = 3a [corresponding to Fig. [TJb)], but 
are still somehow localized, rather at the central part of the GPW for d = 4a [corresponding to 
Fig. Hie)]- These center-localized modes resemble light gradons [corresponding to Figs. Etc) and 
UXf)] at the left tail while look like heavy gradons [corresponding to Figs.[TJa) and[TJd)] in the 
right front, therefore are called "light-heavy gradons." 

Equation © can be mapped into an equivalent chain of graded coupled harmonic os- 



cillators with additional on-site potentials. 8 - In this case, the vibrating mass M n = [1 + 
2e 2 (a; n )] 2 /3Au;pe2(a; n ), and the strength of the additional harmonic spring U n = M n u? n — 2K , 
where K = (a/d ) 3 is a force constant between adjacent masses, A = 2 (—1) for longitudi- 
nal (transverse) case. 9 We then define two characteristic frequencies o; cl (n) = ^JU n /M n and 
w c2 {n) = a/ (Un + AK )/M n . It is easy to notice that cj c i(n) — > u n and u) C 2(n) — > oo n when 
d — > oo, which are as expected because the coupling between the particles vanishes for in- 
finitely large separation. We previously shown that when cj c i(1) < uj c2 (N), one get extended 
modes between uj L (d , c) < lo < uj H (d ,c), where uj L (d , c) = ct> c i(l) and co>h(^o, c) = uo c2 (N). 9 
However, if u> c i(l) > tu c2 (N), there is no extended modes but light-heavy gradons when 
uj L (d ,c) < uo < co> H (d , c), where u L (d , c) = uj c2 (N) and c^h^o, c ) — ^ci(l)- Furthermore, 
there exists a critical point of d = d c when co> cl (l) = to c2 (N), showing only one light-heavy 
gradon which, however, is across the whole system like an extended mode. These arguments work 
well in view of the fact that both the band boundaries and the transition frequencies [cu L (d , c) and 
uJn(do, c)] agree with the numerical data: the two transition frequencies are marked by the vertical 
dashed lines in Figs.|2tb)-|2tf) where they meet with the singularities in the DOS curve (thick dark 
lines). 

As there are various plasmon modes sustained by GPWs, we construct a phase diagram as 
Fig. which shows not only the case of c = 1.0 [Fig.|3tb)] but also the case of c = 0.5 [Fig.|3ta)]. 
Nevertheless, let us focus on Fig. |3jb) for the discussion. Figure |3^b) contains four curves of 
cg> c i(1), u c i(N), cg> c i(1), and to c2 (N) as functions of do- As we have mentioned, there indeed 
exists a point d = d c (vertical dashed line) when cu c i(l) and u; C 2(iV) cross. Also we notice that 
cg> c i(1) = lu c2 (1) — ► loi and u c i(N) = uj c2 (N) — > ujn as d — > oo, which are clearly marked by 
the two horizontal dashed lines of cvi = 0.378^ and u N = 0.333cj p , respectively. In consistent 
with the previous results, we now discuss the four shaded regions partitioned by the four curves in 
Fig-Etb). Specifically, (1) extended mode is possible only when o; cl (l) < uj < uj c2 (N), and when 
do < d c = 3.47a which corresponds to the left panels in Figs, [l] and El (2) light-heavy gradons 
emerge when uj c2 (N) < uj < u; c i(l) for d > d c , corresponding to the right panels in Figs, [l] 
and 12 (3) the lower black region indicates a phase of heavy gradons which have relatively low 
frequencies; (4) the upper dotted region indicates a phase of light gradons which have relatively 
high frequencies. Similarly, the same phase regions exist in Fig. 12a), where a reduced c = 0.5 
results in an increased d c = 4.27a, but the discussions on the four phases through (1)— (4) are still 
applicable. In fact, the critical d c is determined by oj c i(1) = CG> L (d ,c) = uJu(do,c) = to c2 (N) 
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which yields 
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V[l + 2e 2 (x 1 )] 2 ^[l + 2e 2 (x N )] 2 



The results from Eq. © are shown in Fig. |4[d) for both e e = 3.0 and e e = 1.0. This is a guideline 
for determination of extended modes. 15 For instance, only the region below the solid line (shaded 
region) supports extended modes, for the case of ei = 3.0. Note that in this region there exist 
also light and heavy gradons, depending on the frequency, whereas the region above the solid line 
indicates appearance of localized modes only, orderly in the form of heavy gradon, light-heavy 
gradon, and light gradon as frequency increases. 

Up to now, we have already interpreted the mode transitions in view of the equivalent coupled 
harmonic oscillators of graded masses and on-site potentials. From another perspective, we further 
examine the relationship between the site Mie resonance of isolated nanoparticle, the resonant 
band of an infinite PW, and the resonant band of an infinite GPW with infinitesimal gradient. In 
Fig. 12 we have already plotted the corresponding results for the chains in homogeneous host of 
e 2 = eg (dash-dotted lines) and e 2 = q + c (dashed lines), i.e., results for homogeneous PWs. It is 
known that infinite periodic PWs have resonant bands around respective Mie resonance cu n as 4 - 5 



where k is the wave number of the plasmon wave, r the attenuation coefficient. Here the nearest- 
neighboring electrodynamic coupling strength 7^ = jQ n K , where 7 0n is host-dependent. 4 The 
second-order correction due to dissipations is typically cosh(rdo) ~ 1.001 for silver, 4 which is 
negligible. Let us break the GPW into a large number of infinite segments of PWs, each of which 
approximated by homogeneous host e 2 (x n ). For each of these segments we have Eq. ©, which 
means a series of bands located at different Mie resonance u n . Overlapping of these bands can 
be used as guideline to determine mode types, as sketched in Fig. Eta). In Fig. Etb), increased 
d defies the overlapping of these bands, therefore extended modes can not show up. For fixed 
do = 4a, decreasing the gradient coefficient to c = 0.5 will overlap all these bands again, resulting 
in extended modes [see Fig. Etc)]. In this way, the critical interparticle distance d c is determined 
by the condition that the lower band edge of Eq. © for e 2 = e 2 (xi) equals to the upper band edge 
for e 2 = e 2 (x N ), i.e., uj\ — 2A7 2 = lu^ + 2A7^ which simplifies to 



J 1 = — 2A7^ cos(kd ) cosh(rd ) 



(4) 




(5) 
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where 7 0n more explicitly depends on the resonant frequency cu n , the optical effective electron 
mass m*, and the magnitude of the oscillating charge Q. For example, one can write 7 2 n = 
Qe/m* n e 2 {x n )r where e is the electron charge. Note that cu 2 = Qe/rriQ,— we therefore get 

2 3a^e 2 (x n ) = 
l0n l + 2e 2 (x n ) ~ M n X W 

and 

, [l + 2e 2 (x n )] 2 

Here itlq denotes the optical effective electron mass of bulk metal in vacuum. While a free electron 
mass m = 9.1095 x 1CT 28 g, typically = 8.7 x 1CT 28 g for Ag!^ Equation © does not 
account for size and quantum effects in the nanoparticle, however, it is a scale relation which 
indicates that the optical effective electron mass increases from left to the right in the GPWs. This 
is consistent with our terminologies of light, heavy, and light-heavy gradons. Also it is easy to 
show that Eq. © is exactly the same as Eq. ©. 

In conclusion, we have proposed GPW structures which sustain interestingly localized coupled 
plasmon modes that may have potential applications in plasmonics, such as optical switcher and 
multiplexer. The underlying localization mechanism and our analysis equally apply to a wide 
spectrum of analogous problems, which may have important ramifications for understanding ex- 
citations with transitional spectra in many condensed matter systems. 
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Figure Captions 



FIG. 1. (Color online) Real part of the induced dipole moment for longitudinal excitation at 
% = N/2 in a GPW for d = 3a (left panels) and d = Aa (right panels). Damping parameter 
T = O.OOlcjp through (a) - (f). Real frequencies are indicated in each panel. 

FIG. 2. (Color online) Dispersion relation, density of states D(u), and inverse participation ratio 
(IPR) of GPW for d = 3a (left panels) and d = Aa (right panels). Dashed-dotted and dashed 
lines represent the corresponding results for PW in homogeneous host of e 2 = 3.0 and e 2 = 4.0, 
respectively, (a), (d) Dispersion relations for near field coupling, (b), (e) D(uo) versus coupled 
plasmon mode frequency. The thick black lines represent results of GPW with nearest-neighboring 
coupling, (c), (f) IPR versus frequency. 

FIG. 3. (Color online) Phase diagram showing the variation of the four characteristic frequencies 
as functions of d for (a) c = 0.5 and (b) c = 1.0, in the case of = 3.0. The critical interparticle 
distance d c is represented by the vertical dashed lines, whereas the two horizontal lines represent 
uoi and uj n , respectively. 

FIG. 4. (Color online) (a)-(c) Schematics of "band construction" in GPWs. (d) Critical inter- 
particle distance d c as a function of c for both ee = 3.0 (solid line) and e e = 1.0 (dashed line). 
The horizontal dashed line represents d = 2a, which is the geometric limit. Extended modes are 
possible only in the region below these curves, e.g., shaded region in the case of e e = 3.0. 
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FIG. 2: /Xiao, Yakubo, and Yu 
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FIG. 3: /Xiao, Yakubo, and Yu 



(a) d = 3a; c = 1 (b) d = 4a; c = 1 (c) d = 4a; c = 0.5 
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FIG. 4: /Xiao, Yakubo, and Yu 
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